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Abstract 

00 ' 

' The phonon-like excitations of (anti)-instanton liquid (II) due to adiabatic variations 

Q\ ■ of vacuum wave functions are studied in this paper. The kinetic energy term is found 

and proper effective Lagrangian for such excitations is evaluated. The properties of their 
Qh| spectrum, while the corresponding masses are defined by Aqcd with prevailing chromo- 

■ electric component, are investigated based on the phenomenology of QCD vacuum already 

developed. 

(N : 

PACS: 11.15 Kc, 12.38-t, 12.38-Aw 

oo : 

Today's great splash of interest in instanton-based models was provoked by recent attempts to 
Q\ • advance our understanding of QCD behaviour at very high densities [[[]]. It was argued there 

that instantons are able to generate a gap quite suitable to revive the mechanism of the 'colour 
superconductor' in QCD for quantitative utility. In fact, the instantons and, in particular, 
the model of (anti)-instanton liquid correctly seize many nonperturbative phenomena and vac- 
uum features as chiral symmetry breaking, gluon condensate, topological susceptibility [[J>0- 
However, pushing forward the detailed comprehension of QCD principal issues, new questions 
appear that can be addressed to the instanton liquid model for its own sake. We are focused 
on one of those in this paper. 

It is usually supposed the corresponding functional integral in this approach is saturated 
by quasi-classical configurations close to the exact solutions of the Yang-Mills equations (the 
Euclidean solutions called the (anti)-instantons) and the wave function of vacuum being homo- 
geneous in metric space is properly reproduced by averaging over their collective coordinates. 
In the II liquid approach one takes the superposition ansatz of the pseudo-particle (PP) fields 
as a relevant approximation to the 'genuine' vacuum configuration 

N 
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Here A^x;^) denotes the field of a singled (anti)-instanton in singular gauge with 4N C (for 
SU(N C ) group) coordinates 7 = (p, z, fl), of the size p with the coordinate of its centre z and 
Q as its colour orientation. Then, 



Al(x- 1 ) = -n a %^^— , y = x-z, (2) 
9 y 2 y 2 + p 2 
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where 77 is the t'Hooft symbol M. For anti-instanton 77 — 77 (making the choice of the singular 
gauge allows us to sum up the solution preserving the asymptotic behaviour). For simplicity 
we shall not introduce the different symbols for instanton and anti-instanton, and then in the 
superposition of Eq.(P N becomes the total number of PP in the system 4- volume V with the 
density n = N/V. 

The 'liquid' is introduced by the distribution function over instanton sizes and in one-loop 
approximation it is taken as the following 

Hip) = p -B^Nc e -P(p)-^ . ( 3 ) 

As known, the respective parameters are defined by the QCD renormalization scale Aqcd 
together with the constant £ characterizing effective infrared repulsion of the PP's in the model 
of Ref. 0, for example, 

with the Gell-Mann-Low beta function 

P(j>) = - In C Nc - b ln(Ap) , A = = 0.92A P . y . , (5) 
and constant Cn c depending on the regularization scheme, in particular, 

4.6 exp(-1.37iV c ) 
N ^ ^(N c -l)\(N c -2)\ • 

The parameters (3 and (3 are fixed at the characteristic scale p, i.e. 

(3 = f3(p), P = P-\nC Nc . 

The distribution p(p) has obvious physical meaning, namely, the quantity d 4 x dp p(p) is pro- 
portional to the probability to find an instanton of size p in some point of a volume element 
d^x. At small p the behaviour of distribution function is stipulated by the quantum-mechanical 
uncertainty principle preventing a solution being compressed at a point (radiative correction). 
At large p the constraint comes from repulsive interaction between the PP's which is amplified 
with (anti)-instanton size increasing. 

In what follows, we study the excitations of II liquid induced by the adiabatic dilata- 
tional deformations of instanton solutions. We show the corresponding excitations display a 
particle-like ('phonon') behaviour, developing the mass spectrum settled by Aq C d and with the 
chromoelectric component prevailing. It looks quite plausible that incorporating light quark 
condensate would allow us to describe the light hadrons on the correct mass scale, and to 
provide more insight into the confinement mechanism. 

The guiding idea of selecting deformation originates from transparent observation. The 
deformations measured in the units of action (here q, p are the generalized coordinate 
and momentum) have only physical meaning. However, the instantons are characterized by 
'static' coordinates 7 and, therefore, need to appoint the momenta conjugated. It looks quite 
natural for the variable p to introduce p = dp/dx^ as a conjugated momentum. Then we 
should include additional degrees of freedom induced by the variation of instanton size, and 
under sufficiently slow variation of its size the additional terms occur in strength gluon tensor 

Gt = G a + g« . (6) 



2 



The first term corresponds to the contribution generated by the instanton profile 



V ~ 9 {y 2 + P 2 ) 2 l2^° M p y 2 w y 2 J ' 

and in adiabatic approximation the corrections have the form 

dA a 4 y v p 

9a ~ — o P = ~Vaiv , 9 ; nT^ Pi 9ij = > 9iA = ~9a j J ; = 1) 2, 3 . 
dp g {y 2 + p 2 ) 2 

Here we restrict ourselves with the superposition ansatz (|IJ) and neglect the terms of order 0(p 2 ). 
The adiabatic constraint o" . <C Cr" means the variation of instanton size is much smaller than 
characteristic transformation scale of the PP field, p <C 0(1). Then calculating the corrections 
for the action, it is reasonable to take out p beyond the integral and one instanton contribution 
to the action turns out 



5= AG^^ + Cp + ^p 2 , (7) 



1 />_ r /2 _ 87T^ , ri A , K »M ,2 

4i - 9 2 

where p should be taken as the mean rate of slow solution deformation at characteristic instanton 
life-time ~ p. For simplicity, one may take it in the centre of instanton p(z). The constant C 
coming from the interference term should be equal to zero as a result of space homogeneity and 
for the 'kinematical' K-term we have Q 

K s.t.(p) 1 2 3tT 2 / N 



\ J d*x g\ 
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The overt p dependence of k is lacking because of the scale invariance. It arises with the 
renormalization of the coupling constant. Dealing with definition ([I]) we have considered only 
the corrections induced by the variation of strength tensors in line with anzatz but not 
those resulting from a possible variation of fields (0). Bearing in mind the form of potentials 
in regular (r.g.) and singular (s.g.) gauges 



K = -Va^u d u \ia(y 2 + p 2 ), r.g. 



(9) 

s.g. 
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A l = -^Va^ d u hi (l + ^ 
we find the adiabatic corrections = A IM + £is follows: 

a l = ~Van4 -3-7—2 P, r -S- ( 10 ) 

9 y 2 + p 2 

as for the profiles (0), the substitution 77 — > — 77 brings about the transition from regular 
gauge to a singular one. Consequently, it becomes possible to calculate the correction for the 
chromoelectric field (in the singular gauge, for clarity) 



<(P0 = & = ~ g jfj^ Sai , (11) 



In a regular gauge the result is the same. 
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(p. means the potential part here), and the correction for the chromomagnetic field is 



Chromomagnetic field contributes to the action as 

- j d'x KK = -r p\ (is) 

and the contribution of chromoelectric field ( e" = e"(p.) + ef(s.t.) ) is given in the following 
form 

I J W = l ^fp 2 - (14) 
Finally, summing up both contributions to the kinematical K-term results in 

327T 2 . . 

« = — • 15 
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The comparison of Eq. (|13|) and Eq. ([L4]) allows us to conclude that the chromoelectric fields are 
influenced mainly by dilatational deformations. 

The formulae obtained characterize the deformed instanton, and one can see the adiabatic 
variation of the PP solution does not change the action considerably. A small contribution of 
order 0(p 2 ) only is added. Turning to the II liquid, let us rewrite the PP ensemble formulae in 
the compact form of a density functional (as in statistical mechanics). Assuming the number 
of PP in an ensemble is still appropriate to consider them separately, one may rewrite the 
superposition Eq. (P in the form of double summing as 

M x ) = H E M^hn) > (is) 

i=l j=l 

Here AN(pi) is the PP number of the size p in the interval (pi, pi + Ap), K is the number 
of partitions within (pi n , p/ n ) and A IM (x]i, r yj) is the (anti)-instanton solution for PP with the 
calibrated size of p e (pi, pi + Ap) where 7 = (z, ft) stands for the coordinate of centre and 
colour orientation of solution. By definition 2 i=1 AN(pi) = N. The distribution function is 
taken as 

AN(p) 1 

<P) = ~^ v , (17) 

with normalizing Y^i=i n {Pi) Ap V = N. Assuming the limit Ap — > existing one has 
V J dp n(p) = N. As is known, such an ensemble is described by the action averaged over 
the PP coordinates in metric and colour spaces. Then it is easy to obtain (for noninteracting 
system) 

N " diz i jo o r ,a f , / x f 87r2 <p) ^ 



{S ( 0)) = jj| dZi^ s =J d *xJdp n(p) |^ + 



5 We omit the integration over the instanton sizes because we have introduced the collective coordinate 
of density. In principle, the generating functional of the theory should be rewritten in terms of functional 
integration over D[n(p)], however, we use here the distribution function (^) found formerly within one-instanton 
approximation. 
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where the integration over d A z is carried out in the volume V occupied by this system, dQ is 
the measure in colour space with unit normalization (we omit the argument p in the rate p 
hereafter, hoping it is not misleading). This form of (S) reflects the homogeneous property of 
vacuum wave function assumed, together with taking into account possible adiabatic variation 
of the solution profile in various space points. According to Eq.([Tj|) the deformations at a 
constant rate additionally favour non-interacting PP's. 

The situation becomes crucially different if the interaction is switched on. The pair inter- 
action (and the loop corrections) is expected to be of main importance for our purposes. It is 
convenient to set the intensity of gluon field G^ v generated by any pair of PP's as a sum of 
field intensity generated by each separate PP and a term stipulated by nonlinearity of strength 
tensors || 

G>„(a + b) = Gau(a) + G^(b) + AG fll/ (a, b) , 

(19) 

AG p (a,6) = -i{[a^b v \ + [& M ,a„]} , 

(in short G(l + 2) = G 1 + G 2 + G 12 , G^ = G a ^T a and T a are the group generators). The 
cumulative field of the ensemble is given by summing up the integral flTBp and integrate over 
all interacting pairs of PP's 

(s int ) = Jd A x J ^dn, J ^dn 2 ~ x; ^{g 2 12 + 2g,g 2 + 2g,g 12 + 2g 2 g 12 } . 

Colour averaging preserves the first term only and with the explicit form of solution (Q) we 
obtain for the partial contribution of every pair || 

l^ 2N _87r 2 N c f 4 rdS, f d*z 2 \ly\yj - ( Vl y 2 f] p\p\ 

4 l 12) gZN?-lJ dX J V 1 v yKvl + plYvM + pl) 2 ' [ ] 



where yi — x — Zi , pi = pi(zi) , i = 1,2. If the PP sizes do not vary we have for the partial 
contribution 

/Q v 87T 2 2 2 



In the adiabatic approximation the integral (^) is estimated to be (the contact interaction) 

1 . „ . 8tt 2 £ 2 f d 4 z 



Indeed, the adiabatic condition makes it possible to rescale the integration variable — = d (-) + 



&z w d (z) . Besides, it is easy to see that the integral is dominated by the local vicinity of 
diagonal — « 
following form 



diagonal ^ ~ ^ and the distribution function n(p) reduces the interaction term Si n t to the 
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(Sint) = \~^£ 2 J diz J dpl I dp2 Pi U (P^ P 2 U ^ ' ( 21 ) 

remember now pi (i = 1,2) are the functions of z. 

As known, the loop corrections with pre-exponential factor modify the contribution, (|T 
which is estimated in our approach as (in terms of (3 = as defined by Eq. (|)) 

= [d 4 z [ dp n(p) ( ^ p 2 + (3{p) + 51n(Ap) - ln/3 27Vc J . (22) 
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Finally, one has for the averaged action 

(S) = (S®) + J d 4 z I dp n{p) p 2 J dp 2 n(p 2 ) p\ . (23) 

With the distribution function @ for stationary p we easily reproduce the result (§) taking a 
conventional quantity of averaged squared instanton size p 2 = J dp p 2 n(p)/ J dp n(p) . Thus, 
the averaged action is modified to give an effective potential in the following form 

(S) = Jd'zJ dp n{p) p 2 + Ueff{p )\ , 

(24) 

Uefjip) = M + 51n(Ap) - ln/3 2 ^ + v P = . 

p 2 

This result is also used as a suitable approximation herein for the PP's with deformations if it 
includes an additional integration over volume in the definition of p averaged, i.e. 

— Jd 4 zjdp p 2 n{p) 

P Jd^zjdpn(p) ■ { ] 

Everything considered so far concerns the Euclidean action only. In order to get a guess 
for the Minkowski space we notice that in every temporary slice (£4 = const ) we have a 
dynamical system (with a natural transition from p 2 in Euclidean space to — p 2 in Minkowski 
one) where, as easily seen, appreciable benefit in action could occur. In fact, it takes place if the 
rates of solution deformations are self-consistent in every temporary slice and if the oscillator 
fluctuations of excitations are permissible (just the case of interacting potential above). We 
imply no fluctuations of instanton solution itself (their life-time is of instanton size order). Let's 
emphasize that it is worth speaking about the coherent behaviour of vacuum wave function 
when the probabilities of appearance and disappearance of instantons with different rate of 
size deformation are concordant. Besides, the ensemble projection upon the three-dimensional 
space is described by the Lagrangian density for every group of the PP's of calibrated sizes (see 
Eq.flTB]) for partitioning at the initial time) 

C=-\*{p) P 2 + U eff (p) . (26) 

The characteristic scale of excitations is determined at the point of potential minimum 
dUef d t p — = and is derived from the equation pi = ^ p 2 . Then d ^/ ^ \ Pc — =j- Fixing the 
kinetic energy term (p~5| ) on the characteristic scale p (similar to the calculation of potential) 
one has for the frequency of the proper system fluctuations 



m 2 = Jl = = V _ (27) 

k p 2 8n 2 p 2 p(p) p> 



We have only analyzed the deformations in the temporary direction. Those in the spatial 
directions could be estimated by drawing the same arguments. Thus, the expression for the k- 
term keeps the form above obtained with the only change of rates for the appropriate gradients 
of function p(x) done, i.e. the substitution p{t) — > should be performed for such a 

'crumpled' instanton. Then the frequency of proper fluctuations might be interpreted as the 
mass term and excitations occur having a phonon-like nature 

C= l - k( P ) i-p 2 + VpVp] + U eff (p) , (28) 
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(cross-terms ~ p p equal to zero identically), f] The parameters p and (3{p) determined (by 
maximizing the partition function of // liquid with respect to N ||) in a self-consistent way 
take the following values pA « 0.32, f3(p) ~ 18.1, n A 4 « 0.73, [] (iV c = 3) therefore, for mass 
term we have m ~ 1.37A. 

It is a remarkable fact that instanton liquid model gives the effective Lagrangian, incor- 
porating low energy excitations which are naturally interpreted as particles. Their quantum 
numbers are the 'glueball' ones and the mass scale developed is of light hadron mass order. 
Let's emphasize once more that adiabatic assumption leads to fully consistent picture of II 
liquid. The model itself regulates the most suitable regime of dilatation deformations resulting 
in the mass gap generation (|27|). Apparently, an intriguing guess is to associate the lightest 
hadrons with these phonon excitations discovered, since the preliminary evaluations of their 
mass spectrum with the quark condensate included look quite encouraging. The concept of 
confining potential origin for light quarks in this approach seems simply irrelevant because of 
stable phonon nature. 
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6 It is interesting to remark the centre of instanton solution may not be shifted unlike the dilatational mode, 
since the relevant deformations make k singular. 

7 Taken the constant Cn c from Ref. (Nucl. Phys.) as Cn c ~ t'a - P i )7( Af 6 ^2) ! we obtain the following 
values p A w 0.37, (3(p) « 17.5, n A 4 ss 0.44, and for mass gap we have m « 1.21 A. 
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